Probing the Kondo Lattice Model with Alkaline Earth Atoms 
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We study transport properties of alkaline-earth atoms governed by the Kondo Lattice Hamiltonian 
plus a harmonic confining potential, and suggest simple dynamical probes of several different regimes 
of the phase diagram that can be implemented with current experimental techniques. In particular, 
we show how Kondo physics at strong coupling, low density, and in the heavy fermion phase is 
manifest in the dipole oscillations of the conduction band upon displacement of the trap center. 



To date, most cold atom simulations of condensed mat- 
ter systems have focused on the single-band Bose and 
Fermi Hubbard models, both because they are relatively 
simple to simulate and because they arc believed to cap- 
ture a great deal of important physics. However, there 
are many real materials in which the relevance of both 
internal (spin) and external (orbital) electronic degrees 
of freedom preclude description by the single-band Hub- 
bard model. Recently it has been shown that fermionic 
alkaline-earth atoms (AEAs) have unique properties that 
allow for simulation of Hamiltonians with both spin and 
orbital degrees of freedom [I]. Here we discuss avenues 
thereby opened into optical lattice simulations of the 
Kondo Lattice Model (KLM). As is generally the case, 
the necessity to perform cold atom simulation in a trap 
complicates the analogy with the translationally invari- 
ant KLM. However, in this paper we emphasize how 
a trap can help reveal hallmarks of the KLM, includ- 
ing heavy fermion mass enhancement (through dynamics 
induced by trap displacement), and the Kondo insula- 
tor gap (through formation of a density plateau). The 
proposed experimental probing methods (center of mass 
oscillations [2J [3J H] and shell structure density profiles 
[5J |3]) have been demonstrated to be successful diagnos- 
tic tools in alkali atoms and therefore we expect that 
our analysis will have direct applicability in near future 
experiments done with alkaline-earth atoms. 

In its standard form, the KLM consists of a band of 
conduction electrons interacting via a contact Heisenberg 
exchange with a lattice of immobile spins. We focus on 
the case of antiferromagnetic (AF) exchange, relevant to 
so-called heavy fermion materials, known for radically en- 
hanced quasiparticle masses |7J. Simulation of the KLM 
Hamiltonian with AEAs was described in [T], and here 
we only summarize the basic idea. 

The 1 5'o (g) and 3 Po (e) clock states of an AEA can 
be trapped independently in two different optical lattice 
potentials whose periodicities could be engineered to be 
the same [5]. We can therefore consider a Mott insulator 
of e atoms (immobile spins) trapped in a deep optical 
lattice and coexisting with mobile g atoms (conduction 
electrons) trapped in a shallow lattice with the same pe- 
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FIG. 1: (Color Online) [(a)] Schematic zero temperature phase 
diagram for the 1-D KLM [9] ITO] , FM is a ferromagnetic 
phase, and PM is a paramagnetic phase closely related to 
heavy fermions in higher dimensions. Schematic of the ground 
state at strong coupling (b) and for one g atom (c). 
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termined by 4 s-wave scattering lengths a e 
for the states |ee), \gg) and ^ (\eg) ± \ge)) respectively. 
We choose the e atoms to be localized because they would 
otherwise suffer lossy collisions. Collisions between g and 
e atoms on the other hand are expected to be mostly elas- 
tic p] . The independence of the scattering lengths on the 
nuclear spin state guarantees that there will not be spin 
changing collisions, and so we are justified to consider an 
ensemble with only two nuclear spin states a = ± (the 
electron spin in the KLM). 

If the g atoms have negligible interactions with each 
other (true to very good approximation in 171 Yb), and 
the strong repulsion between e atoms is taken into ac- 
count by a unit filling constraint, the low energy Hamil- 
tonian contains only two interaction parameters U t 
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J d 3 rWe (r) Wg (r) (w a is the lowest Wannier orbital 
for the lattice containing the a atoms). Defining V ex = 
(U^g ~ U~ g ) /2, dropping constant terms, and including 
a harmonic trap of curvature f2, the Hamiltonian reduces 
to [I]: 
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FIG. 2: (Color Online) |(a)| In the translationally invariant MFT (fi = 0), hybridization causes the crossing between the 
conduction and localized bands (black dotted lines) to be avoided, and opens a gap Ah- Fillings of n g < 1 correspond to 



filling only the lower band, and mass enhancement is due to placement of the Fermi Surface in the flattened band edges, (b) 
(N g — 25, q — 40, v = 8) The hybridization gap Ah induces a density plateau in the trap. Lines are from self-consistent MFT 
and open shapes from LDA. (c) (N g = 5, q ~ 235) Enhancement of the g atom COM oscillation period (t/to) reveals the mass 
enhancement (m/mo ~ t 2 /tq ! with mo and to the values at v = 0). Circles are from MFT dynamics, curve is a guide to the 
eye. Inset: COM motion for v = 3 (red line) and v = (black dotted line). 
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c laa creates an 



atom at site i in elec- 
{e,g} and (nuclear) spin state a, 
and J g is the tunneling energy for 



In the above 
tronic state a 

the g atoms. The dimensionless ratios v — —2V ex /J g 
and q = AJ g /fl, together with the number of g atoms 
N g , characterize the different parameter regimes of the 
model. It is important to note that both the sign and 
magnitude of V ex will depend on the atomic isotope used, 
and can be adjusted by offsetting one lattice from the 
other (to decrease the overlap between wanier orbitals) 
P]. Therefore, in principle both AF and FM exchange 
are relevant. Nevertheless, in this paper we exclusively 
consider the AF case (v > 0). 

In this regime the phase diagram of the translationally 
invariant KLM in 1-D [Fig. l(a)| has been relatively 
well established, and can be drawn consistently from 
a variety of numerical studies and several exact results 
[3]. At strong coupling ferromagnetism prevails, but the 
weak coupling limit is paramagnetic (PM). The bound- 
ary n g — 1 is insulating, having spin and charge gaps for 
arbitrarily small non-zero coupling. To our knowledge 
the 1-D model is not realized in condensed matter sys- 
tems, but could be explored with AEAs in a 3-D optical 
lattice if both the e and g lattices were made deep in two 
of the dimensions (an array of 1-D tubes). 

We begin our analysis in the PM phase, which is closely 
related to heavy fermion behavior in higher dimensions 
[S]. The mass enhancement can be understood quali- 
tatively through a hybridization mean field decoupling 
[llj in which the g atoms near the Fermi surface gain 
large weight in the localized band [Fig. 2(a) . While the 



mean field theory (MFT) does not capture the Luttingcr 
liquid nature of the PM phase at low energies, we be- 
lieve it nonetheless provides a reasonable guide to the 
phenomena discussed here, and effects beyond MFT are 



left for future study. Moreover, in work to be presented 
elsewhere, the calculations to follow have been extended 
to a 2D geometry where MFT is more reliable, with no 
qualitative change to the results. 

The MFT can be obtained by a (non unique) decou- 
pling of the interaction term in TLk, leading to: 
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In Eq. (2 1 we have denned Vi = \ J2 C 
where the expectation value is taken in the slater deter- 
minant of the (N e + N g )/2 lowest energy single particle 
states (the 1/2 accounts for spin degeneracy). We have 
also introduced chemical potentials /ij to enforce the lo- 
cal constraints (n; e ) = 1- This decoupling is paramag- 
netic, and therefore cannot capture any magnetism, but 
it does describe the tendencies towards singlet formation 
at strong coupling. In addition, it turns out to be the 
exact N — > oo solution of the SU(N) generalization of 
the KLM (which can be implemented with AEAs having 
nuclear spin I = (N — 1) /2) [12j [13j . Beause Hm ft is 
quadratic it can be diagonalized, but it is necessary to 
choose the Vi self consistently. 

In the translationally invariant problem it is custom- 
ary to assume Vi — V and //j = fi, in which case analytic 
progress is possible. With the trap we retain the site de- 
pendent % and Hi, and self consistent solutions must be 
obtained numerically. The procedure involves an initial 
guess for the Vi based on Local Density Approximation 
(LDA): we treat the trap as a site dependent chemical po- 
tential, and infer the energy on each site from a transla- 
tionally invariant problem. LDA results are obtained by 
minimization of the energy thus obtained, while obeying 
a constraint on the total particle number. We then solve 
for fii that satisfy the local constraints [TJ], diagonalize 
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Ti-MFT, and calculate the Vi. Iterating this procedure we 
arrive at a self consistent solution. 

From the MFT ground states we can easily compute 
), which give us density profiles in the trap. For 
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N g or f2 sufficiently large these show plateaus [Fig 2(b) 
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similar to what is observed for the repulsive Hubbard 
Model, although here they reflect the gap of a Kondo 
insulator, not a Mott insulator. The Kondo insulator 
can be understood from the hybridized band picture [Fig 
2(a)] . Unit filling of g atoms correspongs to completely 
filling the lower band, and there is a charge gap of Ah- 
LDA considerations then imply that ft (jf — jf) = Ah 
2(b)| . Exact results for the v = oo KLM give 
- jf) = 3\V ex \; in this limit A H tends to 2 \V ex \, 
so the MFT underestimates the plateau size. For the 
bosonic Hubbard model, where the relevant gap is the on- 
site interaction U, such plateau structures have already 
been imaged via microwave spectroscopy 0[6]. We there- 
fore expect that for large v the plateau can be observable 
experimentally. 

At lower fillings, such that the plateau does not form, 
we are everywhere in the heavy fermion metallic state. 
Under these conditions we consider an experiment where 
the trap center is suddenly displaced, causing dipole os- 
cillations of the g atom center of mass (COM). These 
type of experiments have been implemented in alkali 
atoms to study 1-D and 3-D transport of interacting 
bosons and fermions [3J [31 H] and used to probe differ- 
ent quantum many-body regimes in these systems. We 
calculate these dynamics self consistently, starting with 
the MFT ground states and shifting the Hamiltonian. 

If a \a = E 4 ( u i4g<r + vl q 4ea) create the eigenstates of 
H-mft then the following set of discrete Schrodinger 
equations governs the evolution of the u q (t) and v q (t) 
after displacing the trap by S lattice sites: 

-*«? = * + U ? +1 ) - n (i - sf u\ - v ex Vivi 

= -n{i-5) 2 vl - mvl -V ex %u\. (3) 

When integrating these equations the V are updated self 
consistently, and the /ii are evolved in time by ensuring 
that h e i — (the first time derivative has no dependence 
on the Hi ) . Such time dependence of the \ii preserves the 
one e atom per site constraint, and its necessity has a 
simple origin: Hmft breaks the local U(l) symmetry of 
Hk associated with conservation of the e atom density. 

The g atom COM oscillations ensuing from displace- 
ment of the trap by one lattice site have been obtained at 
several different values of v for fixed q. We find a strong 
enhancement of the oscillation period t (and hence of 
the quasiparticle mass m ~ t 2 ) for decreasing \V ex \ [Fig. 
2(c)| . Once v ~ 1, the comparably fast non-interacting 
oscillations emerge on top of the slow oscillations of the 
heavy quasiparticles, converging to the noninteracting re- 
sult as v — > - this cuts off the apparent trend toward 
diverging r in Fig. 2(c)| 
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FIG. 3: (Color Online) (a) Non interacting dynamics of 16 g 



atoms after displacement of the trap by 5 lattice sites. 
Same as|(a)| but now with v = oo (calculated using Hoo 
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We now turn to the FM part of the phase diagram, 
which exists for all fillings at sufficiently strong coupling. 
When v — oo the ground state is formed by pairing each 
g atom into a spin singlet with one e atom [Fig. [jjb)]. To 
first order in J g the singlets become mobile, and there is 
an exact mapping of the unpaired e atoms to the fermions 
of a U = oo Hubbard model (the singlets are the holes) 
[15] . Because nearest neighbor hopping cannot exchange 
up and down spins, we can think of the fermions as spin- 
less [16] . perform a particle hole transformation (now the 
singlets arc the spinlcss fermions), and thereby arrive at 
a simple Hamiltonian describing the g atoms: 
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The reduction of the hopping energy is the result of pro- 
jecting out the high energy triplet states. To highlight 
this strong coupling behavior, we again consider dynam- 
ics ensuing from a sudden displacement of the trap cen- 
ter. If the system is strongly interacting, Eq. Q avails a 
simple treatment of these dynamics based on the non- 
interacting solutions in Ref. [T7]. There the authors 
found that for q ^> 1 and N g < 4-^/2 J g /Q the dynamics 
involved delocalized, free space harmonic oscillator like 
states with level spacing to* = Q^/q. At v = oo we ef- 
fectively have N g — > 2N g (because the fermions become 
spinlcss) and J g — > J g /2, therefore the inequality can 
be violated at strong coupling even when satisfied for 
the non-interacting system; localized states become pop- 
ulated, and transport is strongly inhibited (see Fig. [3]). 

Another limit which is well understood in the FM 
phase of the translationally invariant model is N g = 1. 
Sigrist et al. [18 proved that the ground state of the 
KLM with L sites and one conduction electron has total 
spin S = \ (L — 1) even in the absence of translational 
symmetry, and they described the excitations for the ho- 
mogenous case as bound states between the g atom and 
a flipped spin in the deep lattice (spin polarons, see Fig. 
[IJb)). For a weak trap (q ^> 1) and at sufficiently small 
coupling (V ex -C to*), we characterize the polaron spec- 
trum to lowest order in degenerate perturbation theory 
and find that one eigenvalue separates from the rest by 
a gap of approximately \V ex \ [Fig. 4(a)| . As seen in Fig. 
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[4(b)] this energy scale is manifest in the COM oscillations 
of a single g atom as strong modulation of the oscillation 
amplitude with periodicity r w rfrn; as verified by dy- 
namics calculated from the exact eigenstates. 
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q = 800, showing a gap of size 



(Color Online) (a) Low energ y po laron spectrum for 
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Polaron COM (black dotted) oscillating after trap displace- 
ment by 2 lattice sites. For comparison we plot cos(V ex t) 
(blue solid line), showing that the energy gap determines the 
time scale of amplitude modulation. The overall decay is a 
finite size effect, and eventually revives. 

For the finite system under consideration, we expect 
this behavior will persist for N g > 1 whenever N g is odd. 
This can be seen by noting that a single unpaired g atom 
at the Fermi level gains energy at first order in pertur- 
bation theory when coupling to the e atoms is turned 
on, whereas the doubly occupied levels below it do so 
only at second order. It is worth noting that the con- 
dition V ex <C lu* is equivalent to demanding the pertur- 
bation stay smaller than the finite size gap. For fixed 
lattice depth and peak g atom density, the gap scaling 
is w* ~ with R the Thomas Fermi radius of the 

g atom cloud (R cx ^ N g J g /oj*). This means that the 
demonstrated modulations will be washed out with in- 
creasing \V ex \ or with increasing g atom number, and are 
manifestly a finite size effect. 

We now consider the feasibility of generating and ob- 
serving these dynamics in an experiment. Throughout 
the paper we assume a unit filled Mott insultor of e atoms 
coexisting with various fillings of g atoms at the center 
of a trap. To realize this situation in experiment, one 
could first ramp up a deep optical lattice for the g atoms 
such that they exhibit a Mott insulator shell structure. 
By taking advantage of the energy shift between dou- 
ble and single occupied sites, it is then possible to selec- 
tively transfer atoms into the e state, in such a way that 
sites with two gg atoms become eg sites and single occu- 
pied g sites become single occupied e sites. A subsequent 
adiabatic reduction of the lattice depth for the g atoms 
achieves the desired configuration. It may also be helpful 
to confine the g atoms more tightly than the e atoms (to 
ensure that they do not sample the wings of the e atom 
Mott insulator), which can be achieved by blue detuning 
the deep lattice. In most of our calculations we have used 
small trap displacements to simplify the numerics, and in 
a real experiment they will inevitably be larger. Obser- 
vation of the dipole oscillations in 1-D with amplitude of 



< 8 lattice sites has precedent, and was carried out by 
mapping the center of mass position of the atomic cloud 
to momentum space[21[3]. Moreover, a recent proposal 
[TO] suggests that dynamics could be characterized from 
a single non destructive measurement if the atoms are 
coupled to an unpumped cavity field. The small atom 
numbers necessary for observation of the modulations at 
small | | also has precedent, with ~ 18 atoms per tube 
having been realized in an array of ID lattices [50] • We 
emphasize that all physics discussed in the paper involves 
energy scales on the order of V ex , which makes tempera- 
ture demands less constraining than for proposals involv- 
ing superexchange or RKKY like physics (v ss 1 gives 
a Kondo temperature on the order of \V ex \ /ks, so this 
statement applies even to the heavy fermion behavior). 

The results presented demonstrate that dipole oscil- 
lations of the g atom COM effectively probe a variety 
of KLM phenomena. The emphasis has been on a ID 
system, which is of relevance to cold atom experiments, 
primarily because current theoretical understanding of 
the phase diagram is strongest here. However, the heavy 
fermion behavior generalizes to experiments in 2 and 3 di- 
mensions. Though it has not been discussed, we point out 
that an optical lattice experiment — especially in D > 1 — 
is a natural setting for probing the size of the Fermi sur- 
face in the heavy fermion state, and could corroborate 
evidence for a large Fermi surface observed in condensed 
matter experiments. 
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